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Particle production by slow-changing gravitational fields is usually described using quantum field 
theory in curved spacetime. Calculations require a definition of the vacuum state, which can be 
given using the adiabatic (WKB) approximation. I investigate the best attainable precision of 
the resulting approximate definition of the particle number. The standard WKB ansatz yields a 
divergent asymptotic series in the adiabatic parameter. I derive a novel formula for the optimal 
number of terms in that series and demonstrate that the error of the optimally truncated WKB 
series is exponentially small. This precision is still insufficient to describe particle production from 
vacuum, which is typically also exponentially small. An adequately precise approximation can be 
found by improving the WKB ansatz through perturbation theory. I show quantitatively that the 
fundamentally unavoidable imprecision in the definition of particle number in a time-dependent 
background is equal to the particle production expected to occur during that epoch. The results are 
illustrated by analytic and numerical examples. 



I. INTRODUCTION 

Particle production by gravity in a slowly expanding 
universe can be described usingquantunr field theory in 
curved spacetime (QFTCS) HtH S S H . A well- 
known feature of QFTCS is the absence of an absolute 
definition of vacuum and particles for quantum fields in 
arbitrary curved backgrounds (see e.g. 0). The num- 
ber of particles detected by an observer depends on the 
observer's motion, and there are no preferred observers 
in a general spacetime. If the universe is sufficiently spa- 
tially flat and expands sufficiently slowly so that the four- 
curvature scale is much larger than the wavelength of a 
field mode, there is a natural (if approximate) definition 
of the vacuum state for that mode: the adiabatic vac- 
uum with respect to a given fiducial time to. This is the 
vacuum state seen by "approximately inertial" observers 
at t = to- However, an adiabatic vacuum state defined 
at t = to is generally an excited state with respect to the 
vacuum defined at another time t\ ^ to. Such field exci- 
tations are interpreted as particles produced by gravity. 
Since there is no single physically preferred vacuum state, 
we speak of an apparent particle production. The "real" 
particle content of a given quantum state of the field can- 
not be unambiguously established without referring to a 
particular physical experiment where the vacuum state 
is observed or prepared. 

For instance, the most straightforward definition of the 
vacuum state (the "instantaneous diagonalization" of the 
Hamiltonian) yields an infinite apparent particle produc- 
tion in some generic cases |j , while the adiabatic vacuum 
exhibits finite particle densities in the same cases. It has 
been also proposed ^(| that the vacuum state should 
be chosen to minimize the apparent particle number ob- 
served at a time to- This prescription is physically rea- 
sonable and yields a quantum state close to an adiabatic 
vacuum, but the resulting state will in general depend on 
the choice of to- A natural and unambiguous definition 
of the vacuum state is available only in regimes where 
gravity becomes negligible. So the concept of particles 



can be used in curved spacetimes only in an approximate 
sense which becomes more precise in slow- changing, al- 
most flat geometries, and also for high-energy particles. 
The central theme of this paper is to explicitly analyze 
the precision of this approximation. 

The prescription of the adiabatic vacuum is based on 
the WKB approximation and has been particularly useful 
in the context of QFTCS (see e.g. [ul GH E3). It is 
well known that the WKB approximation is applicable 
to equations such as 

^L+^(t)x = 0, (1) 

where uj(t) > is a time-dependent frequency function. 
The approximate solutions can be found in the form of an 
asymptotic series in the adiabatic parameter T , where 
T is the characteristic variation timescale of the function 
Lu(t), which is assumed to be slow-changing 13]. Us- 
ing the n-th order WKB approximation, one can define 
adiabatic vacuum states of order n The difference 

between n-th order adiabatic vacua defined at different 
fiducial times is characterized by apparent particle occu- 
pation numbers, which are of order this is then 
the imprecision in the resulting definition of particle num- 
bers. A brief overview of the calculation of cosmological 
particle production using the adiabatic vacuum prescrip- 
tion is given in Sec. I11BI 

It is known that the WKB approximation involves a 
divergent asymptotic series 1 and thus cannot be used be- 
yond a certain order n max - Therefore, particle occupation 
numbers computed with respect to an adiabatic vacuum 
are defined only with a certain fundamentally limited ac- 
curacy. If one uses the adiabatic vacuum of the optimal 



I shall show in Sec. HTDl that the WKB series generally diverges. 
Although this statement appears to be common knowledge, I was 
unable to find a derivation in the literature. A closely related 
result is th e B orel summability of asymptotic series for adiabatic 
invariants ll-ll : see also Ref. Il5l . 
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order n max , one obtains particle occupation numbers up 
to an uncertainty of order T ~ nnLax , and this accuracy 
cannot be improved any further. 

Since typical particle numbers produced in vacuum by 
smooth geometries are exponentially small |2J, namely 
of order exp(— uT), it is a priori unclear whether even 
the best attainable precision ~ r~ 1_ ™ ma * is adequate for 
the description of particle production. The exponentially 
small particle numbers can be calculated if one applies 
perturbation theory techniques to the WKB ansatz; var- 
ious such techniques are outlined in Sec. 11111 I shall sum- 
marily refer to these techniques as the perturbatively im- 
proved WKB. 

Although the WKB expansion, being a power series in 
T _1 , necessarily misses any exponentially small terms, 
the quantity T _,lmax can be of order exp(— u>T) if n max is 
sufficiently large, say of order T. The main result of the 
present work is an explicit estimate of the optimal order 
"max and of the resulting optimal precision of the WKB 
series. The current literature does not appear to offer 
such direct estimates. It is difficult to analyze the WKB 
series since there is no closed-form expression for the n-th 
term of that series. To circumvent this difficulty, I use 
a particular perturbatively improved WKB technique — 
the so-called Bremmer series (see Sec. 1111 J?|) . I show in 
Sec. llVl that n max — 0(u>T) and the error of the optimally 
truncated series is exponentially small. More precisely, 
for Eq. Q with analytic functions u>(t), the optimal order 
"max of the WKB approximation at t — to is found as 



u(t)dt 



(2) 



where the complex numbers ti, i = 1,2,... are all the 
zeros and the poles of u>(t) in the complex t plane, and 
the integrals are performed along the paths from to to tj 
that give the smallest value to the above integrals. The 
error of the optimally truncated WKB series is of order 



max 



exp(-2n max ) . 



(3) 



I also show that the smallest attainable uncertainty in the 
definition of particle numbers during a time-dependent 
epoch is of the order of the particle production expected 
to occur within that epoch. These issues are discussed 
in Sec. II V CI I illustrate these estimates by analytic and 
numerical examples in Sec. 



II. ADIABATIC APPROXIMATION IN QFTCS 



A. The WKB approximation 



where e is a formal parameter (we shall set e = 1 at the 
end of all calculations). The frequency u)(t) is assumed 
to be a sufficiently slow- varying function of time, so that 
the adiabaticity condition 

du> 



dt 



(5) 



holds for all values of time t to be considered. We shall 
additionally assume throughout this paper that iv 2 (t) > 
for all relevant t, and that u>(t) is an analytic function. 
The well-known WKB ansatz is 

xvnsB{t) = ^==exp (±- f uj{t')dt'\ (6) 



The error of this approximation is of order e 2 if ui(t) is 
a sufficiently well-behaved function pH ll7T| . One can 
generalize the ansatz (0 to an n-th order approximation 

see e.g. El El El), 



(n) 
C WKB 



(*) 



c 



exp ± 



W n (t')dt' 



(7) 



where the auxiliary function W n {t;e) is an approximate 
solution of 
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3 W' 2 
4 W 2 



W 2 



-2 



which is found as a power series in e 

'1 u> 3uj 2 



Y^e 2 * Sk {t)=u-e 2 [~-~ 
^ V 4 uj 2 8 to 3 

k=0 
1 UU^ 



(8) 



(9) 



16 uj 4 



5 cJu> 



13 

32 uj 5 



99 u)lo 2 
32~ 



297 to 3 
128^ 



The approximation a;{^ B (t) is accurate up to error terms 
of order 0(e 2n+1 ), so the series ©, whether convergent 
or not, is an asymptotic expansion at e — > |l3j. The 
consecutive terms of the series @ can be found itera- 
tively, for instance, by expanding the RHS of the relation 



+i 



0(e 



2n+4\ _ 



1 w;i 3 w< 2 

2 W n 4 W 2 



(10) 



in powers of e 2 up to order 0(e 2n+2 ). It is clear that 
W n (t;e) is a rational function of ui(t) and its derivatives 
up to the order 2n. No closed- form expression appears 
to be available for the n-th term S n (t) of the series 

The WKB expansion can be equivalently parametrized 
by the characteristic timescale T of the variation of u>(t). 
After rescaling the time variable by t = Tt, the WKB 
series becomes a power series in T~ x . However, we shall 
not use this parametrization. 



The WKB approximation, also known in the mathe- 
matical literature as the phase integral method and the 
Liouville- Green approximation pH Il7|. applies to equa- 
tions of the form JT|) . rewritten as 



,d x 



+ L) 2 {t)x = 0, 



(4) 



2 See, however, Ref. 12(1 for some attempts to simplify the WKB 
terms by eliminating total derivatives, Ref. Il'll for an approach 
to make high-order WKB calculations numerically more manage- 
able, and Ref. |22| for an approximate resummation using Airy 
functions. 
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B. Quantum fields in expanding universe 

In this section I briefly review the computational pro- 
cedure for determining the (apparent) particle occupa- 
tion numbers in QFTCS, following 7]. To be specific, 
I consider a minimally coupled massive scalar field (f>(x) 
in a Friedmann- Robertson- Walker universe with flat spa- 
tial sections and the line element dr 2 — a 2 (r)(ix 2 , where 
a(r) is the scale factor assumed to be a known function 
of time r. It is convenient to pass to the conformal time 
t = J r a _1 (r)rfr (which is denoted by t for consistency 
with the previous notation) and to rescale the field <j> as 
4>{x) = a~ l x{x). The auxiliary field xi x ) is quantized 
using a mode expansion of the form 



d 3 k 1 



(a k e lk x z; fc (i) + H.c) , (11) 



where Ok are the annihilation operators, Vk(t) are mode 
functions for the wavenumber k, and "iJ.c." denotes the 
Hermitian conjugate terms. The mode functions Vk(t) 
are complex- valued solutions of the equation 



+ [k 2 + 



a" 

| Vk 

a 



v k 



- Vk 



subject to the normalization condition 



Im {v k v* k ' 



1, 



0, (12) 



(13) 



where m is the mass of the field <fi and the prime ' denotes 
derivatives with respect to t. (In a flat spacetime, we 
would have u>f.(t) = const and v k ~ e -i^ k tfe ^ rj ue t 
Eq. (| 1 31) . the creation and annihilation operators satisfy 
the standard commutation relations, [fit, a^,] = 5(k— k'). 
The vacuum state |0) of the field is defined as usual by 



a k |0) = for all k. 



(14) 



Thus defined, the vacuum state depends on the choice 
of the mode functions Vk(t). Because of the freedom to 
multiply each mode by a constant phase factor, solutions 
of Eqs. (|12|) - (|13fl may be effectively parametrized by the 
(complex- valued) ratio v' k /vk at a fixed time t = t$. Dif- 
ferent choices of this ratio yield mode functions v k (i) de- 
scribing different vacuum states. 

We shall now focus on the behavior of one field mode 
with a fixed k and hence drop the subscript k. The next 
step is to apply the WKB approximation to Eq. i|12[l ■ The 
n-th order WKB approximation to the mode function is 



vwkbW 



exp 



W n (t')dt'), (15) 



where the factor y/e ensures that the normalization con- 
dition ifTl^l) holds. Using the function «wkb(*), one de- 
fines the adiabatic vacuum of order n at a fiducial time 
t = to by requiring that the mode function v(t) should 
match the WKB expression at t = to, namely 



v'(to) = ^wkb^o) 
v(to) v WKB (t )' 



(16) 



Let us denote the resulting mode function by vo(t). 
(More generally, one may require that the matching in 
Eq. I|16|l should hold only up to terms of order e 2n+2 , i.e. 
up to the precision of the n-th order WKB approximant, 
but we shall not make use of this additional freedom.) 

The adiabatic vacuum prescription can be applied at 
a different fiducial time t = t±, yielding another mode 
function v\ (t) . These two mode functions are related by 
a Bogolyubov transformation 



v {t) = av 1 (t) + Pv* 1 (t). 



(17) 



It is well known that the t — to vacuum appears to have 
the number density \f3\ 2 of particles with respect to the 
t = t\ vacuum. Once the functions vo(t) and vi(t) are 
found, the Bogolyubov coefficient (3 can be computed as 



v'o{t)vi{t) -vojty^t) 
2i 



(18) 



The r.h.s. of Eq. i|18|) is a time-independent Wronskian 
and thus can be evaluated at arbitrary t, say at t = t\. 
Since the values of «i(£i) and i>i(£i) are known from 
Eq. (|16f> after replacing to by t\, it remains to compute 
vo{t\) and v' (ti). The latter computation requires solv- 
ing Eq. I|12|) through the interval [to, ti] with initial condi- 
tions (|16fl . Note that the n-th order WKB approximation 
to vo(t) identically satisfies the condition (|16fl for all t and 
thus cannot be used to determine the particle number at 
t = t\. The mode function i>o(i) must be computed ei- 
ther by a more accurate analytic method or numerically. 
In principle, a higher-order WKB approximant could be 
used if its precision were adequate, but this is not always 
the case, as we shall show in the next subsection. (The 
value of P obtained from a higher-order WKB approxi- 
mant would be of order e 2n + 2 while the correct answer is 
usually exponentially small.) 



C. Super-adiabatic regimes 

Let us consider the case where oj(t) sa const except for 
a finite time interval, for instance uj = ujq for t < to and 
lu = uji for t > t\ (the values to and t\ may be finite 
or infinite). A possible such function uj(t) is plotted in 
Fig. ^ I n th & t case, the vacuum states are naturally 
and uniquely defined for t < to (the "in" vacuum) and for 
t > t\ (the "out" vacuum). We shall refer to this situation 
as an "in-out" transition, and to the regimes t < to and 
t > t\ as super-adiabatic regimes. 

More precisely, a super-adiabatic regime at t = to 
means that the adiabaticity condition (jSJ) becomes an 
equality, 



lim 



0. 



(19) 



and that analogous conditions hold for higher derivatives 
of uj(t). In other words, all derivatives of uj(t) vanish in 
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Figure 1: A frequency function u>(t) with two super-adiabatic 
regimes at t < to and t>t\. 



special choices of tu(t). Below we shall obtain more pre- 
cise estimates of the growth of terms in that series, and 
presently we demonstrate the generic divergence of the 
WKB series using qualitative arguments. 

Let us consider a frequency function co(t) that allows 
an "in-out" transition between t = t$ and t = t%, with 
super-adiabatic regimes at t < to and t > t\. Suppose 
that the series (0 converges for some e = Eq and for 
all t within the relevant range [to? *i] - Then this series 
will converge absolutely for smaller |e| < Eq, yielding a 
well-defined function 



a super-adiabatic regime. The WKB series is truncated, 
so that we have 



W n (t) = 



u>q for t <to, 
u>i for t > t\, 



(20) 



because W n (t) is a local function of uj(t) and its deriva- 
tives at a point t, as seen from Eq. @. Thus, within 
one super-adiabatic regime, the adiabatic vacuum states 
of all orders coincide and yield a natural definition of the 
vacuum state and an unambiguous notion of particles. 
For instance, the result of an "in-out" transition is a well- 
defined set of particle numbers (the particle density of the 
"in"- vacuum with respect to the "out"- vacuum) . Physi- 
cally this means that gravity becomes unimportant in a 
super-adiabatic regime, and all inertial detectors exactly 
agree on the particle content of any quantum state of the 
field. Outside of a super-adiabatic regime, adiabatic vac- 
uum states are still well-defined, but there can be only 
an approximate agreement between adiabatic vacua of 
different orders or defined at different fiducial times. 

If the quantum field is in the "in"-vacuum state, the 
particle numbers observed at t > t\ can be unambigu- 
ously predicted using Eq. I|18l) . However, if we used the 
WKB approximation to compute the function vo(t) at 
t = t±, we would find vo(t) oc e IWlt / £ , the Wronskian in- 
volved in Eq. (|18|l would vanish, and we would obtain an 
incorrect result (3 = 0. It is well known that the particle 
number is generically nonzero, except for certain special 
cases when the particle production exactly vanishes. 

We conclude that the WKB ansatz approximates the 
mode function vo(t) insufficiently accurately for cal- 
culations of particle numbers. I shall outline some 
known methods of improving the WKB approximation 
in Sec. Ell 



D. Divergence of the WKB series 

If the function oj(t) is analytic (or at least C°°) in t, we 
may attempt to sum the series (0 by evaluating the limit 
n — ► oo at a fixed t. However, the series (0 generally will 
not converge as n — > oo. Although the initial terms may 
decrease, eventually after a large enough n the terms S n 
will grow without bound. Thus the series © can be in- 
terpreted only as an asymptotic series, except for certain 



fc=0 



s k (t). 



(21) 



By construction, the function Woo(t; e) is analytic in e at 
e = 0. We shall now show that Woo is an exact solution 
of Eq. ©. Since the partial sums W n were obtained by 
expanding the recurrence relation 11(J|) to a finite order 
in e, we may substitute Woo instead of W n in the r.h.s. of 
Eq. 10 and find that 



W n + 0(e 2n+2 ) 



IW" 3W£ 



2 Woo 4W£ 
holds for all n. Therefore, the relation 



(22) 



Woo = 



2 Woo 



3 W' 2 

_ vr oo 

4 W 2 



(23) 



holds as an identity between power series in e after both 
sides are fully expanded. However, the r.h.s. of Eq. IL'-'SI) 
can be also viewed as a power series in e in which only 
the square root has been expanded, namely 



Wo 



2lu 



IW" 3W r 



2 W a 



4 W 2 

^ vv oo 



(24) 



Since the l.h.s. of Eq. I)23|l is an absolutely convergent 
series, and the terms of such a series can be reordered to 
form the r.h.s. of Eq. (|24ll . both series specify the same 
analytic function of e. Hence, the series on the r.h.s. of 
Eq. (|24|) converges and thus the function Woo(t;s) is an 
exact solution of Eq. (|23|l . 

It follows that the "infinite-order" WKB mode function 



Voo(t) 



VWo7{t) 



exp 



Woo (t')dt' 



(25) 



is an exact solution of Eq. Then, from the defini- 
tion of the adiabatic vacuum state, we see that the mode 
function Voo(t) represents the adiabatic vacuum of ev- 
ery order and for all fiducial times t at once. In this 
uniquely identified adiabatic vacuum state, the particle 
number is exactly zero at all times and thus there is no 
(apparent) particle production at any time. This out- 
come contradicts the fact that particle production in vac- 
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uum is present in the generic case. Therefore in general 
the "infinite-order" WKB solution ^(i) cannot exist. 3 

Another way to arrive at the same conclusion is to con- 
sider the particle number |/3(i;e) 2 | in the second super- 
adiabatic regime (t > ti) as a function of e. We can 
express |/3(t;e)| through the function Woo{t\e); for ex- 
ample, using a zeroth-order adiabatic vacuum, we find 



\P(t-e)Y 



e 



(u- W c 



\2uj 



W 

* * PC 

2W„ 



(26) 

i oi2 



Note that the rapidly oscillating phase is absent from |/ 
and that Woo(t] e = 0) = w(t) for all t. Thus \j3(t; e)\ 2 is 
nonsingular (and equal to zero) at e = 0. On the other 
hand, the particle number is known to decay faster than 
any power of e if the function u>(t) is analytic or smooth 
(of class C°°). Thus the particle number in the second 
super-adiabatic regime, |/3(ti; e)| 2 , is an analytic function 
of e which decays faster than any power of e at e = 0. 
Such a function must be identically equal to zero. It fol- 
lows that the particle number \f3(ti;e)\ 2 would be iden- 
tically zero for all \e\ < £ max , given the convergence as- 
sumption 1)21(1. We know, however, that generic choices 
of Lj{t) involving an "in-out" transition will exhibit non- 
vanishing particle production 5 even for very small e. We 
conclude that the WKB series cannot converge for all t 
within the interval [£o,ii], except for special choices of 
uj{t) where the particle production is absent. (Note that 
the WKB series converges to u>(t) for t within a super- 
adiabatic regime since all derivatives of u>(t) vanish there, 
but nevertheless it does not yield a solution x{t) valid for 
all other t. We shall later show that in general the WKB 
series cannot converge even for some special values of t 
outside of super-adiabatic regimes.) 

Finally, let us illustrate the growth of terms Sk in the 
series @ by performing an estimate of only the highest- 
order time derivatives of co(t) entering the expressions 
Sk- This estimate is merely qualitative because in fact 
the contributions of lower-order derivatives cannot be ne- 
glected; we shall obtain more precise results in Sec. IIVI 



Rewriting the recurrence relation (|1U|> as 




(27) 



we find after some algebra that the term S n contains the 
highest-order derivative of u>(t) always in the following 
combination, 



S n (t) 



(-1)' 



1 



n,,,2n 



A n u) 



(lower-order), (28) 



where we have suppressed terms containing lower-order 
derivatives of uj. Now we need to estimate the growth of 
derivatives of For a generic analytic function f(z) 

having some poles in the complex z plane, the growth of 
derivatives at z — zq can be estimated as 



d n f(z) 



dz r ' 



(Zi - zo 



.n+l ' 



(29) 



where Z\ is the pole of f(z) nearest to z = zq (see Ap- 
pendix ^ for a derivation of this formula). In our case, 
poles of l/V'*' correspond to zeros of u>(t), i.e. turning 
points in the complex t plane. Let us suppose that the 
nearest zero of u(t) is at t = t%. Under the assump- 
tion that the value of S n is of order of the first term in 
Eq. QgJ, we find 



\Sn(t)\ 



2u(t) 



(2n)! 



\t - 1 



i2n+l • 



(30) 



It is now straightforward to see that the terms S n may 
decrease for small n but eventually begin to grow after 
n* ~ e _1 w(i) \t — ti\. The term S nt is of order 



2d, 



(2ra.)! 



1 



\t-ti 



i2n»+l 



exp(— 2n*) . (31) 



The above values are roughly correct order-of-magnitude 
estimates for the optimal number of terms in the WKB 
series and for the optimal precision, as we show below. 



III. PERTURB ATIVE IMPROVEMENT OF 
WKB 



3 Of course, an exact solution Wooif) of Eq. 1231 exists, but it can- 
not be exactly represented by the WKB series J5J because such 
Woo (t) must be quickly oscillating and cannot be a slow-changing 
function of t, as implicitly assumed by the expansion 1^1 . 

4 Kulsrud llgl has shown that the change in the adiabatic invariant 
is of order e p+2 if u>(t) has p continuous derivatives. One can 
also demonstrate an exponential-like behavior |/3| oc exp(— C/e) 
for analytic u>(t), under some technical conditions 23]. See also 
Ref. |24| for a derivation of an exact invariant which coincides 
with the adiabatic invariant in the leading order in e. 

5 Exceptions are provided by cases where all t erms S„(t) of the 
WKB series vanish after some n (see e.g. l25l . chapter 2, for an 
explicit example), and also by functions uj(t) corresponding to 
Schrodinger equations with exactly zero scattering |2(I . 



The precision of the WKB approximation can be im- 
proved by applying a time-dependent perturbation the- 
ory to the WKB ansatz. Equivalent and closely related 
perturbative methods were used e.g. in |3.lll| for particle 
production calculations and in the mathematical litera- 
ture [l6L fill for precision estimates of the WKB ansatz. 
A related technique called "quasilinearization" [23 was 
applied in e.g. |28l |29| to obtain improved approxima- 
tions to bound states of the Schrodinger equation (al- 
though the convergence of the method was investigated 
only numerically). I shall now give a brief overview of 
these methods. 
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A. Methods based on the Riccati equation 

The simplest version of the perturbative improvement 
technique is intended to determine a small correction to 
an approximate positive- frequency solution of Eq. , 



x+ (t) = C exp 



i i uj(t')dt' 
>o 



in the form 



x(t) = x+it) exp 



y(t')dt' 



(32) 



(33) 



Thus one introduces a new (complex-valued) dependent 
variable y(t) via 



X 
X 



iuj(t). 



(34) 



The equation for y(t) is easily derived and is a Riccati 
equation, 



y — 2iujy = id) — y 



(35) 



By assumption, x(t) differs little from x+(t), and so one 
expects that the value of y(t) is small and that the y 2 
term in Eq. (|35|) can be treated as a perturbation. Dis- 
regarding the y 2 term and using the natural initial con- 
dition y(0) = 0, one finds the approximate solution 



?/(!)(*) 



i I ui(t') exp 
/o 



li I cj(t")df 



dt'. 



(36) 



This first approximation y\ (t) is already sufficient to com- 
pute the leading contribution to the exponentially small 
particle occupation number in the context of an "in-out" 
transition. The resulting (approximate) value of (3 is 
found after some algebra as 



(37) 



For analytic functions w(i), the value of /3(i) (ti ) is typi- 
cally exponentially small due to rapid oscillations of the 
integrand in Eq. i|36|) . 

A variation of this perturbative improvement tech- 
nique was developed and used in Ref. [3(|. There the 
new dependent variable for Eq. was defined by 



X — IUJX 

The function Q{t) satisfies the equation 



C-2^C = --(i-C 2 ) 



(38) 



(39) 



and can be interpreted as the instantaneous squeezing 
parameter describing the squeezed state of one mode of 
the quantum field at time t relative to the instantaneous 



vacuum state at that time. The particle number is found 
from 



W)\ = 



icwr 



i - icwr 



(40) 



The leading-order solution Qi) of Eq. (|39|l with the initial 
condition £(0) = is obtained by disregarding £ 2 , 



C(i)(*) = - 



exp 



2i / uj{t")dt" 



(41) 



Further approximations C(n) (t) are found from the recur- 
rence relation 



C(n+i)(*) - _ 



jdt' 



2u> 



i-Cf n) (*') 



exp 



2i / ui(t")dt" 



(42) 

It was shown in Ref. |3pJ that the sequence C(2) ; ■•• 
converges to the solution of Eq. as long as | C( n ) (t) \ < 
1. An advantage of using the variable instead of y(t) 
is that the values of £ are always bounded as long as 
w(t) 0, which helps in the analysis. (It is straightfor- 
ward to verify that Eq. (|13|) together with Im u> = 0, ui =^ 
yield the bound \(\ < 1.) 

The method of "quasilinearization" [23 uses the ansatz 



x (2 ) (*) = C exp 



(l/(2) (t) +2/(i) (t) -iw) dt' 



which leads to the easily solvable equation 



y ( 2) - 2 (iu - 2/(1)) 2/(2) 



(43) 



(44) 



where the term quadratic in y^ 2 ) nas been neglected and 
2/(i) is given by Eq. I|3t)|l . Further corrections 2/(3), 2/(4) 5 ■■• 
and the corresponding solutions X( 3 ) , x^ , ... are deter- 
mined in the same manner. Analytic and numerical stud- 
ies |28| show that the accuracy of the solution X(„)(i) im- 
proves quadratically with n, the error being of order e 2 . 
However, precise conditions for the convergence of this 
method do not seem to have been investigated. 

The methods outlined so far are conceptually simple 
but lead to nonlinear equations which complicates their 
analysis. We shall therefore use an equivalent but some- 
what more long-winded perturbative technique based on 
a system of two linear equations. The advantage will be 
that we shall obtain approximations to the solution x(t) 
in the form of a series (called the Bremmer series). 



The Bremmer series 



In this section I mostly follow Refs. |31|,|32j. The usual 
WKB approximation to solutions of Eq. Q is 



x(t) w AX + (t) + BI_(t), 
where A, B are constants and 
1 



X ± (t) 



exp 



1 



u(t')dt' 



lo 



(45) 



(46) 
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are the two WKB branches (here, to is an arbitrary initial 
point). To improve this approximation, one looks for 
solutions x(t) in the form 



x(t)=p(t)X + (t)+q(t)X_(t), 



(47) 



where the coefficients p{t), q(t) are now time-dependent. 
By introducing two unknown functions p(t), q(t) instead 
of one unknown x(t), we have added a degree of freedom 
which is canceled by imposing the additional relation 

x(t) = l -tu(t) [-p(t)X+(t) + q(f)X- (t)] . (48) 

This relation signifies that the derivative of x(t) can be 
computed from Eqs. (|46(l - (|47|l by formally treating p(t), 
q(t), and Lo(t) as constants. It is then straightforward to 
derive the following simple equations for p(t) and q(t), 



(49) 



In the context of QFTCS, one is usually looking for 
small corrections to the positive-frequency WKB solu- 
tion x{t) — X+{t). Therefore we assume the initial con- 
ditions p(to) = 1, q(to) = 0, and rewrite Eqs. (|49|l as the 
following integral equations, 

P(t) = l+o/ -l^rft, q(t) = - / --±pdt. (50) 



2 Jt w X 4 



2 7 t0 ^X_ J 



These equations can be solved iteratively, starting from 
the initial approximation p(t) = 1, q(t) = 0. The result 
is conveniently written as a series involving an auxiliary 
sequence u n (t), 

00 00 

p(t) = l + J2 u 2n{t), q{t) = J2 u 2n-l{t), (51) 
n— 1 n— 1 

where the functions u n (t), n>0, are defined recursively 
starting from uo(t) = 1 as 

M 1 

U2n(t) = - / — -T^U2n-ldt, H > 1, (52) 

(53) 

The resulting infinite series for the solution x(t), 

oo 

x(t) = X + + Y,(u 2n ^ 1 X^+u 2 nX + ), (54) 

n=l 

is called the Brcmmer series. A physical motivation be- 
hind this derivation is given in [23j and references therein. 

C. Convergence of the Bremmer series 



, , i r uo x+ 

VQn+i{t) = - ——u 2n dt, n>0. 

Z J t U A - 



t max may be infinite) under the rather weak restriction 6 

dt < oo. (55) 

The convergence can be demonstrated by the following 
argument pH l32l l34| . The sequence u n (t) is majorized 
by the auxiliary sequence U n (t) defined by 



r 


Co 


'to 


LO 



U n+ i{t) = 





Co 


f 


2ui 



U n dt, U (t) = 1; 



(56) 



namely < U n (t) for all n, t. Convergence of the 

majorizing series Y^ =0 U n (t) is therefore sufficient for 
the absolute convergence of the Bremmer series. In turn, 
the series z~2^=o U n (t) = U(t) can be summed explicitly 
by deriving a differential equation for U(t), 



LO 

2^ 



U(t) 
The solution is 

U(t) = cxp 



U(t), U(t ) = 1. 



to 



2uj 



dt 



(57) 



(58) 



Thus the Bremmer series converges as long as the integral 
in Eq. I)58|l is finite, which yields the condition l|55(l . The 
convergence is uniform in t because an upper bound 





LO 


r 

'to 


2lo 



dt < M 



entails 



\u n (t)\ < U n (t) < -M n , t <t<t n 



(59) 



(60) 



and so the number of terms n may be chosen in advance 
to guarantee a desired precision for all t < £ max . Namely, 
it is straightforwardly seen that a relative precision e~ p 
will be achieved by computing n = 0(P/ \nP) terms. 



IV. PRECISION OF THE WKB 
APPROXIMATION 

A. The adiabatic expansion 

The WKB series (jHJ entails an adiabatic expansion of 
the WKB solution 



^wkbW = 



1 



exp 



W n (t')dt' 



to 



(61) 



The Bremmer series converges to the exact solution 
x(t) absolutely and uniformly for to < t < i max (where 



6 The condition 1551 would be violated, for instance, in the pres- 
ence of parametric resonance with t max = oo. 
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in powers of the adiabatic parameter 



^WKB(i) = X + (t) 



is f ( u> 3 lO 2 \ , 
4 L„ \lu 2 2uj 3 



e 1 ( u 3 u? 

T V^ 5 ~ 2^ 



(62) 



where X + (i) is defined by Eq. (|46|l . The scries in brackets 
above contains terms produced by expanding the denom- 
inator l/y/Wn as well as by expanding the exponential. 
Therefore we expect that the series in Eq. (|62() is again 
a divergent asymptotic series. We shall now analyze this 
series to determine the optimal number of terms and the 
best attainable precision. 

The main idea of this analysis is to represent the ex- 
act solution x(t) by the (convergent) Bremmer series and 
to compare the latter with the WKB series l|62|) . How- 
ever, the structure of the Bremmer series differs from that 
of the WKB series in two aspects. Firstly, the Brem- 
mer series involves both positive- and negative-frequency 
branches while the WKB ansatz H62|) contains only the 
positive-frequency branch. Secondly, the Bremmer series 
is not a power series in e, but instead it involves oscillat- 
ing integrals containing e -1 under exponentials. Thus we 
need to obtain an asymptotic expansion of the Bremmer 
series ansatz l|47|) in powers of e, of the form 

oo oo 

p(t)X+ + q{t)X_ =X+Y, A k e k + X_ B ^ k , (63) 



fc=0 



k=0 



where A k and B k are some time-dependent coefficients. 
The positive-frequency part (X + ^ fc Ak£ k ) of the above 
expansion will coincide with the expansion (|62() because 
the asymptotic expansion in power series is unique [35| . 
Note that the split between positive- and negative- 
frequency branches in the expansion in the r.h.s. of 
Eq. H63|) is well-defined [ljj, but nevertheless ^Z fe Ak£ k 
is not an expansion of pit), and neither is ^2 k Bk£ k an 
expansion of q(t). We shall now derive an asymptotic ex- 
pansion of the form (|tj3|> from Eqs. (|52fl - l|53l) . It will turn 
out that only the positive-frequency branch survives in 
the asymptotic expansion, i.e. B k = 0, which makes the 
correspondence with the WKB series immediate. 

Let us pass from the time variable t to the dimension- 
less "phase" variable 9, 



9{t) 



u(t')dt'. 



(64) 



This is a well-defined transformation as long as u>{t) re- 
mains real and nonzero in the relevant range of t, and in 
that case ui(t) remains an analytic function of 9. Equa- 
tions 1|52|) - (|53[1 are then rewritten as 



where for convenience we have introduced the dimension- 
less function 



a(6) 



1 duj 
2^~dt 



(67) 



Note that ea(6) -€i I under the adiabaticity condition JSJ. 

Starting from Eqs. I|65|l - H66(l with uq(9) = 1, we shall 
first obtain expansions for U\(&) and u 2 {9), which will 
make further calculations more transparent. To expand 
ui(9) in powers of e, we repeatedly apply integration by 
parts to Eq. I|66|) with n = and find 



ui{0) 



(2if (2i)- 



'=0 

(68) 

where the prime denotes derivatives with respect to 9. 
We further assume that the point 9 = (t = to) is located 
within a super- adiabatic regime where a(9) and all its 
derivatives vanish. Then Eq. I|68|) simplifies to 



ui(9) 



-2i6/e 



2i 



E 



d k a{9) 



d9 k 



k =o (20 

The expansion for 11-2(6) is found from Eq. H65|) as 



(69) 



u 2 (9) 



t x t 



(70) 



note the absence of the quickly oscillating exponential. 
Hence, the initial terms of the expansion l|63l) are 



X + (l + u 2 ) + X- Ul =X A 



is 



~7i 



3tf 
2~ 4 



)dt 



(71) 

which reproduces the initial terms of Eq. (|62|) . The 
negative-frequency branch is absent due to the factor 
exp(— 2i9/e) in u\. It is easy to see that further terms u n 
of the Bremmer series have expansions of the same form: 
namely, the odd-numbered terms u 2 n+i contain the oscil- 
lating factor while the even- numbered terms u 2n do not. 
For convenience, we separate these factors explicitly and 
define the auxiliary functions X n (9) and fi n (9) by 

w 2n -i(0)= (-^) n e- 2te/e \ n (9), »>1; (72) 
u 2n (9) = (--) ix n {6), n>0. (73) 

The functions X n , \j, n are (formal) series in e that satisfy 
the following recurrence relations, 

e\ k d k 



u 2n (9) = f e 2ie '/ e a(6')u2n-i(e'W, n>l, (65) A »+iW = E ( J) ^ K^KW] , n > 0; (74) 

JO k=0 
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u 2n +i(9)= e- 2ie l e a(Q')u 2n {O')d0', n > 0, (66) 



a(9')X n (9')d9' , n > 1; /i = 1. (75) 
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In terms of these functions, the asymptotic expansion of 
the Bremmer series can be expressed as 

oo 

p(6)X+(9) + q(6)X_ (9) = X + (9) A k e k 



k=0 



fe=0 



(76) 



where we have formally defined A = 0. Since the r.h.s. of 
Eq. (0 will coincide with the WKB series (|62[l after A/c 
and /ifc are fully expanded in powers of e, it is sufficient to 
analyze the convergence of these latter expansions. Note 
that the n-th order WKB approximation corresponds to 
retaining the terms of the series (f76*jl up to k = 2n. 



B. Precision of the asymptotic expansion 

The goal of this section is to derive the optimal trunca- 
tion of the asymptotic series l|76ll . We begin by analyzing 
just the first two terms of that expansion, namely 



(u + u 2 )X + + uiX- = A% 



1-^CAx+mi) ■ (77) 



Subsequently we shall examine the convergence of the 
terms \k,^k with k > 1. 

We find from Eqs. I|69|) and l|7U|l that the quantities Ai 
and /ii are expressed by the series 



« (2i) 

,k 



k d Qk ' 



fe =o ( 2z ) Jo 



a{9>) 



,,d k a{9') jQt 



d9' k 



d0' 



(78) 



(79) 



These series are typically divergent because derivatives 
d k a/d9 k of an analytic function a(9) grow as A;! with 
k — > oo. The growth of such derivatives is determined by 
the location of the singularities of a(9) in the complex 9 
plane (see Appendix 1X1 for more details). For instance, if 
the singularity of a(9) nearest to the point 9 is a simple 
pole at 9 1 with residue —c±, then for large k there is an 
asymptotic estimate 



d k a(9) 
d9 k 



c x k\ 



(80) 



In the present case, the function a(9) defined by Eq. i|6T|) 
has simple poles at the locations 9i,9 2 ,--- corresponding 
to zeros t±, t 2 , •■■ of u>(t), namely 



w(t)dt, j = 1,2, 



U>(tj 



0. 



11) 



A pole of u(t) corresponds to 9 = 00 and typically will not 
be the nearest singularity of a (6). If u>(t) has a zero of v- 
th order at t = tj, i.e. ui(t) ~ (t — tj) v : the corresponding 



residue of the simple pole of a(9) will be 



1 



2^+1 



(82) 



Assuming that 9\ is the nearest pole of a(9), it is 
straightforward to find that the terms of the series l|69() 
will start growing after 



k > K(9;e) = 2s- 1 \9 1 -6\= 2e 



uj{t')dt' 



(83) 



If there is any ambiguity in the choice of the complex 
integration contours in Eqs. I|81|) and (|83|) . the contours 
should be such as to yield the smallest value of \9\ — 9\ 
since we require 9\ to be the nearest singularity to 9. 

The best attainable precision 5u\ is of order of the 
smallest retained term, which can be estimated using 

/ 2nke~ k k k , as 



Eq. I|8U|) and Stirling's formula, k\ 



5u\ 



(!) 





d k *a(9) 




d9 k * 




-k, 



(84) 



We shall show below that higher terms k > 3, 
of the Bremmer series have better convergence behavior 
and their error is dominated by that given in Eq. I|84|) . 
Thus the formulae l|83|) - (|84|l are the central result of the 
present paper. The order n of the WKB approximation, 
as defined in Sec. Ill AI is related to the order k of the 
corresponding Bremmer series by 



~-k*(t;e) 



uj{t')dt' 



(85) 



Hence, the optimal order of the WKB approximation and 
the optimal precision are given by Eqs. (J2J and |3"|. where 
we have set e = 1. The prefactor c\ is of order unity and 
will be relatively unimportant for our considerations. 

Let us now analyze the convergence of the series (|7U|I . 
After repeated integration by parts, we have 



(86) 



and the first term dominates for large k. Thus we find 
the following estimate for large k, 



« a(0)a<*- x >(0) « a(9) Cl{k ~ 1) '\ 



(87) 



Hence, the asymptotic series l|70|) starts to diverge after 
k > 1 + fc*(0; e), i.e. one term later than the series for u\, 
and the best attainable precision Su2 is 



Su 2 = (-) a(9) 



d k *a{9) 



d9 k - 



2 V 



(88) 

Note that 6u 2 = (ea/2) 8u\, while ea 1 due to the 
adiabaticity condition, hence 8u% <C 8u±. We find that 
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the precision of the partial asymptotic series Ai + [i\ is 
limited by the error in Ax, which is estimated by Eq. (|84p. 

It remains to analyze the behavior of the higher-order 
terms in the expansion 176( 1 . It follows from Eq. I|74|l that 
the series for A„ is always a sum of several series of the 
form 

E(s) w m - m 

k=Q 

where f(9) is some analytic function expressed through 
a(9) and its derivatives. The functions f(9) will have 
singularities at the same points in the complex 9 plane 
as the function a{9). Thus each series of the form 1(89(1 
will have essentially the same convergence properties as 
the series for \\{9) analyzed above; namely, the series will 
start diverging after fc* terms and have an exponentially 
small error term. However, in the expansion 1)76(1 the 
series X n is multiplied by e". Thus, the series for U2n-i, 
where n > 1, always converges better (i.e. it starts to 
diverge at a later term) than the series ((69(1 for u\, and 
the lowest error in U2 n -i is always smaller than that of U\. 
The same considerations can be seen to apply to the series 
for fx n , n > 1, and hence for the terms Ui n , n > 1, of the 
Bremmer series. Therefore, the optimal truncation of the 
expansion GEJ is the same as the optimal truncation of 
the series Xi(9), as given by Eq. i(83|) . and the lowest error 
is dominated by the error in Ai, as estimated by Eq. (|54"ll . 
This argument concludes the derivation of Eqs. ©-©■ 



where 9\ is the pole of a(9) closest to 9 = 9q in the com- 
plex 9 plane, and the function a{9) is defined by Eq. (|6"7Jl . 
Now consider the total particle number produced during 
the active epoch. As reviewed in Sec. lIIII the leading con- 
tribution to this particle number is exponentially small 
in the adiabatic parameter and can be computed from 
a perturbatively improved WKB approximation. For in- 
stance, using Eqs. l(4"U|) and (|4"T)l yields 

/ + OC 
e~ 2ie ' e a{9)d9. (91) 
-oo 

Since by assumption the function a(0) is analytic, the 
integral involved in Eq. ((91(1 can be transformed into a 
contour integral in the lower half-plane. The latter in- 
tegral contains contributions from all poles 9j of a(9), 
weighted by e~ 2l9j ^ e . The dominant contribution comes 
from the pole which is closest to the real line, say 9\. 
Heuristically, this pole will be also the closest to the ac- 
tive epoch around 9 — 9q. Hence 

/3~exp(-2e- 1 Im(0 o -0 1 )), (92) 

which is of the same order as exp (— 2n max (0 o ))- The pre- 
cision of the WKB approximation away from 9 = 9q will 
be better than exp(— 2n max (6>o)), however, the particle 
production will also decrease. We conclude that the un- 
certainty inherent in the definition of particles during an 
"active epoch" is quantitatively the same as the intensity 
of particle production at that time. 



C. Particle production and the definition of 
vacuum 

We have seen that the WKB series diverges unless 
there is exactly no particle production. Furthermore, we 
shall now show that the ambiguity in the definition of 
vacuum is quantitatively related to the presence of par- 
ticle production. 

Let us consider an "in-out" transition between two 
super-adiabatic epochs and assume that there is nonzero 
particle creation, as there will be in the generic case. 
Intuitively one expects that most of the particles are cre- 
ated during the epoch where gravity is most significant 
(the "active" epoch). However, the adiabatic vacuum dur- 
ing the active epoch is defined only up to the precision 
of the WKB approximation which is fundamentally lim- 
ited. The precision is improved as we move away from 
the active epoch, but so does the expected rate of parti- 
cle production. We shall now derive some estimates that 
illustrate this connection. 

As we have seen, the optimal precision of the WKB 
approximation is of the order exp(— 2n max ), where n max 
depends on the time at which the WKB approximation is 
applied. We shall use the variable 9 defined by Eq. l(6"4"|) 
as the time variable. If the active epoch corresponds to 
values around 9 = 9 , then, according to Eq. l(55|) . 



(90) 



V. EXAMPLES 

In this section we shall investigate the precision of the 
WKB approximation for particular functions uj{t). We 
use the formalism developed in the previous sections and 
set e = 1. 

The first example is 



u>(t) = u>q I 1 + A tanh 



(93) 



where loojA and T are constants. There are two super- 
adiabatic regimes at t — > ±oo. Let us compute the opti- 
mal order n max for the WKB approximation ^(vkb(^o) 
some intermediate time to- The function uj(t) has zeros 
at 



t = t* + iirkT, k £ Z; 



-T 



2 1 



.4 



A 



(94) 



The closest zeros to a real point to will be it* . Then the 
optimal order of the WKB series is found as 



){t)dt 



u; 



t* — t(] 



, „ , cosh 
AT In T 



which grows linearly with to for to - 

"max ~ WO |*0 I (1 ± A) = Lo(t ) | *0 | 



cosh ^ 

(95) 

±00, namely 
to -> ±00. (96) 
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Figure 2: Magnitudes of first 10 terms S n , n — 1, 10, of 
the WKB series for u(t) given by Eq. I|930 and different val- 
ues of to- The values S n are normalized to So and shown in 
logarithmic scale. The curves are labeled 1-8 according to the 
consecutive values of to in Table Crosses indicate the error 
estimates; abscissas of crosses are set to the predicted values 
of n m ax(to)- The crosses are close to the lowest-error terms 
on the corresponding curves, indicating a good agreement be- 
tween the predicted and the numerically obtained values. 



The smallest term of the series is estimated as 



■ exp (-2n max ) 



(97) 



The worst precision is found near to ~ where we have 
(for A < 1) 



: wo 



U + AT In cosh ^ 



TTLUqT 



For the purposes of numerical calculation, we chose 
u! = 3, A = 2, and T = 3/2. The magnitudes of the first 
10 terms Si, Sio of the WKB series © are plotted in 
Fig. |2 together with the error estimates (see Table [I] for 
numerical values). It is clear from the plot that the terms 
S n indeed start to grow after about n = n max , and that 
the error estimate agrees with the numerically obtained 
smallest terms of the series. 

For comparison, let us approximately compute the par- 
ticle occupation number after the "in-out" transition from 
t = — oo to t — +oo. The leading-order contribution to 
the particle number is proportional to the term ui of the 
Bremmer series. It follows from Eq. (|53f) that 



lim m {t) 



at — exp 

2uj 



2i / oj(t')dt' . (99) 



to 


-2.0 


-1.5 


-1.0 


-0.5 


0.0 


0.5 


1.0 


1.5 


^raax 


1.5 


1.2 


1.0 


1.2 


2.0 


3.7 


6.0 


8.4 


s„ m „ 

Sn 


0.06 


0.12 


0.17 


0.08 


0.007 


10~ 4 


HT 6 


6-10~ 9 



Table I: Numerical values of the predicted optimal truncation 
order n max and error estimates 5„ max (see Fig. The pre- 
dicted values of n max are not integer since they are computed 
from Eq. (gSJ. 



The integral can be evaluated as a sum of residues of the 
integrand over its complex poles at 



t = U-mkT, k = 0,1, 2,. 



(100) 



Only the simple poles in the lower half-plane need to be 
accounted for. (The function uj(t) also has poles at t = 
(«7r/2 + ink) T, which do not contribute to the integral.) 
After some algebra, the result is found to be 



lim u\(t) 

t — *+oo 



2sinh [ttw T(1 - A)] ' 



(101) 



Thus the particle number is of the same order as the worst 
precision of the WKB series, ~ exp(— 2n max ), where n max 
is estimated by Eq. I|98|) for A <C 1. This confirms the 
general conclusions of Sec. IIVCI 
The second example is 



Lu(t) = LJ ( 1 + 



(102) 



The four complex roots of u)(t) are 

t fc =Texp^ + ^, k = 0,1,2,3. (103) 

The optimal truncation order is estimated as the smallest 
of the following four numbers, 



mm 

k 







/ uj(t)dt 


= luq min 


J t 


k 



to + 



(tk ~ to) 



5T 4 



(104) 

For t > 0, the value of n max is determined by the closest 
roots ±Texp(i7r/4). For instance, at t = we have 
"max ~ ^oT, while for large t > T we get n max w 
^woto? 1-4 . The rapid growth of the allowed number of 
terms for to — > oo indicates a super-adiabatic regime. 

For the numerical calculation, we chose loq = 1 and 
T = 2. The results are shown in Fig. Eland Table HJ 

The final example involves the function 



w(t) 




1*1 < T; 



(105) 



which exhibits super-adiabatic regimes at t < —T and 
t > T. We shall assume that \A\ < 1. The function uj(t) 
has zeros at 



t k =±TJl 



1 



- In A + 2ink ' 



k € Z, 



(106) 
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Figure 3: Analogous plot as in Fig. for the function ui(t) 
given by Eq. 11U2H . The curves are labeled 1-8 according to 
the consecutive values of to in Table [H] 
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Table II: Numerical values of the predicted optimal truncation 
order and error for Fig. [3] 



which form sequences converging to the essential singu- 
larities at t — ±T (see Fig. . Thus the nearest singu- 
larity to a point to on the real line is effectively the point 
t = T . Then the estimate © gives 



}{t)dt 



In 



lo (T- t ), 0<i <T. (107) 



However, this estimate was derived only for functions 
with simple isolated zeros, while in the present case 
there exist infinitely many zeros in any neighborhood of 
t = ±T. So we do not expect a close agreement with the 
numerically obtained values of n max - 

For the numerical calculation, we chose u>q = 1, T = 
8, and A = 0.1. The results are given in Fig. [5] and 
Table IIIII The computed values show that a few first 
terms of the WKB series can be used near to = T; the 
estimate n max « is invalid in that regime. 



VI. SUMMARY 




Figure 5: Analogous plot as in Fig. |2]for the function ui(t) 
given by Eq. (110511 . The curves are labeled 1-7 according 
to the consecutive values of to in Table 11111 Relatively poor 
agreement between crosses (the estimated values of n max ) and 
the numerically obtained values is due to the presence of es- 
sential singularities in cu(t). 
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In this paper I have reviewed some aspects of the adi- 
abatic approximation and its application to cosmological 
particle creation, and presented new results. It is well 



Table III: Numerical values of the optimal truncation order 
and error for Fig.Q)] 
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known that there exists a fundamental limit on the accu- 
racy of the notion of particles in curved spacetimes, due 
to Heisenberg uncertainty relations Q. A quantitative 
investigation of this accuracy is the main focus of the 
present paper. I showed that the best attainable preci- 
sion in the definition of particles is exponentially small 
and of the same order as the typical particle production 
expected during the same epoch. The conclusion is that 
the ambiguity inherent in the definition of particles is 
precisely due to the possibility of particle production. 

The main technical issue was to obtain an explicit es- 
timate of the highest attainable precision of the WKB 
approximation. I have demonstrated that the WKB ap- 
proximation involves a divergent series and derived a 
novel formula 10 for the optimal truncation order ~i max 
of that series. I also estimated the error of the optimally 
truncated WKB series [Eq. ©]. The error is exponen- 
tially small since typical values of n max will be large if the 
adiabaticity condition JSJ holds. Physically, the value of 
"max is determined by the number of oscillation periods 
during the time of appreciable change in the frequency 
Lu(t). Finally, I have presented analytic and numerical 
examples illustrating the validity of these estimates. 
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and it is well known that this series converges absolutely 
within a circle 1 2; — ^0 1 < Ri, where R\ is the distance 
between the point 9q and the nearest singularity of w(9) 
in the complex 9 plane. Suppose for simplicity that the 
function w(9) has simple poles at 9 = 9j with residues 
— Cj, and that 9\ is the pole nearest to 9q. Then we may 
express w(9) as 



w{9) 



Cl 



(A2) 



where the auxiliary function wi(9) is analytic within a 
larger circle | z — 9 \ < R 2 , R2 > Ri ■ The n-th derivative 
of w(9) is therefore 



w (n \0 o ) 



m ™\ n+1 +v>?\e ). (A3) 
(9i - oq) 

The same reasoning may be applied to the function wi (9) 
and the result is the following asymptotic estimate for the 
growth of derivatives, 



w {n) (9 ) = n\ 



Cl 



C2 



-Oo) 



n+1 



\n+l 



n 

1?! 



— Cl [i + o((iW#i) n )] 



(A4) 



I am indebted to Gerald Dunne, Larry Ford, Stephen 
Fulling, Matthew Parry, and Alex Vilenkin for helpful 
discussions. 



The estimate (|A4|) is straightforwardly generalized to 
the case when w(9) has poles of higher order, e.g. 



w{9) 



Cl 



w 2 (9), 



Appendix A: GROWTH OF DERIVATIVES OF 
ANALYTIC FUNCTIONS 



in which case we have 



Here I derive some formulae used in Sec. lIIDl and llVI 

By definition, an analytic function w(9) can be repre- 
sented by a Taylor series at a regular point 9q, 



71 = 



(Al) 



w {n) (9 Q ) ~ ^ ((n + 1)^ + ci) [1 + O ((R1/R2T)} ■ 

Similar estimates can be easily obtained for the case when 
w(9) has more than one pole at the same distance from 
the initial point 9$. 
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